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ABSTRACT 


ое Osma ais thesis is to investigate the state 
of the art of radar signal design as well as radar signal 
processors and determine the actual trends in modern radar 
design. The use of a digital general purpose radar signal 
MmeeecssOr 15 discussed. The concepts of ambiguity and auto- 
Ewrelation function are investigated in regard to radar 
m olution Capabilities. Тһе concept and analytical devel- 
opment of the DFT/FFT are presented. Quantization noise in 
По са processor and its effects in the improvement 
Ио are analyzed. Optimization techniques for the response 
curve of digital MTI processors using staggered PRF are in- 
vestigated. The SAR concept and analysis as well as techniques 
ENObtaln low correlator rates in the SAR digital processors 


ОТО presented. 
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ПИО ОО СЛОМ 


The processing of signals has been always a necessity 
whenever information was to be transmitted through a channel 
between a source and a destination. It was not until 1948 
e hannon presented the basic theory of information re- 
(што the source entropy with channel capacity and the 
Mmeemapility of error. Techniques of signal processing were 
EMEGdy in use, but the results given by the information 
шивогу clarified the limits that could be achieved. 

Early analog processors with enough bandwidth and reason- 
пе о топа! to noise ratio could process data at extremely 
ПІ Пата rates, but the components' low stability and high 
ШІ апа sometimes low versatility, were responsible for a 
Euxdual substitution of analog by digital processors when 
Nu EDr1ce of digital logic went down, and the digital hard- 
Pewee technology suffered a big jump in the last ten years. 

ЕЕ. тп today's radar signal processing the trend to 
ПВ ЕСИ Стол of analog by digital processors is inevitable. 
Mé theory and boundaries of digital signal processing are 
continually being improved as well as hardware components. 
Ecdgscalmost no longer a problem and versatility is immense. 

Ewan introduction to radar history. In Chapter 
111 signal processing is related with the radar problem through 
the radar range equation and a general digital processor for 


a radar system is presented. In Chapter IV an analysis is 





madero rte present signal processing techniques in radar 
нереста emphasis into digital approaches. In Chapter V 
ШШ ілусетірастоп is made into the recent developments and 
Peoplcms Of twomemayor areas of radar digital signal process- 
Bern Chapter an overall view of the achievements and 


Meenas Of radar Signal processing is explored. 





ткен ГЗТОКҮ AND APPLICATIONS 


The electromagnetism and electromagnetic wave propagation 
leo ries ot Maxwell, later (1886) experimentally verified by 
Hu. contain all the necessary background to understand the 
Demet pies Of radar. The first experiments on detection with 
NAO waves started in 1905. Due to inadequate technology 
Peul ts were very poor. The obtained ranges were less 
meme cnose achieved with optical systems and this was enough 
MmomaascOurage any interest in pursuing the experiments. It 
feemonly in 1922 at the Naval Research Laboratory that a 
wooden ship was detected using a CW radar with the transmitter 
ее receiver as separate units. From then on an increas- 
Memerncerest in radar technology became evident. In 1930 
Ма bistatic radar, the first aircraft was detected, апа 
ШЕ |9252 the detection ranges were already in the 50-mile 
ШЕШІП. Бу that time the pulse techniques were not yet 
Істей, so the information was in the presence or absence 
EM target and not in range information. But in 1935 the 
ИЙЕ ЕП successfully used pulse techniques to measure distance 
Pe targets. Іп 1936 the detection range was already in 
tre M niles Пош ie baste principles were understood. 
Preeiiext Steps were in technology improvements in order to 
ns table elements and more power output to increase the 
ieee coverage. Higher frequencies were used in order to 


work with smaller elements and high gain antennas. 


tO 





The discovery of the magnetron was an important step since 
powers went up by a factor of 100, and wavelengths of 10 cm 
Ex scbtaurned But the greatest incentive for further devel- 
opments in the radar field came from the military necessities 
Or World War II. After the War a stagnation of about five 
years slowed down the rhythm of improvements in the field. 
РІП (Пс Liities the introduction of the high power klystron 
E ucssed not only the power output but the frequency stability 
EN essary for coherent detection as well as coherent MTI; also 
EN low noise receivers were implemented. 

Paoli them On, With more and more perfect technology, more 
tex systems appeared. The introduction of the computer 
EN U-torape and control element was of extreme importance. 
All ballistic systems and very large coverage systems suffered 
Ко impulse. Digital systems improved the implementation 
iS tereable array antennas, synthetic aperture radar, track- 
ing systems, as well as all digital processing techniques. 

With the development of radar, the statistical nature of 
the detection problem gave rise to an increasing interest in 
BIN:tudy of the statistical properties of clutter, radar 
EN section, rain, etc. So, due to a more precise knowledge 
(modeling) of the systems environment, the radar systems 
ed depending on the type of application. 

Шіне, а асстоп to Optimize a given application seldom 
moe USed jin Others. Today the radar design field is wider. 
ШОШ Сагтоп5 Vary from military to civilian to scientific. 
imp rovements in technology and more accurate modeling, 


ШІатгсп іп the field is far from becoming saturated. 
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III. SIGNAL PROCESSING IN RADAR 


A. INTRODUCTION TO THE RADAR RANGE EQUATION 


Only the analysis of the equation 


К = f(n,,n5, ... nj) (5-1) 
Пеште К 15 the distance from the radar antenna to the target, 
and (n4,n5, e Dy) are the radar and environmental parameters, 
ive the necessary knowledge to maximize R. If P. 9125 121016 
meee pulse power radiated and С the gain of the antenna, the 


pususmitted power density is 


Р, © 





Дт К^ 


AE received power density, after reflection in a target 


DE TOss section M will be 


рб о 
t bs (3-2) 


(41R?)* 


(5-2) і5 multiplied by the antenna aperture Аа» the 


received power will be 


Р A 
р. EN AU Ее (3-3) 
AR) 
eg n is the minimum received signal, then the maximum 
distance is 
PECATU 
E Ве (151) 
2 о 
(Ап) аа 


2 





ЕО сегітттіоп оғ по15е figure t Oí a receiver 


© 
= 4 0 1 = 3 - 
5. N; m ( No ) i = input (3-5) 
O2 OU EDU 
where 
N; = КТВ, = input noise (3-6) 


Then, combining (3-4) with (3-5) and (3-6) 
_ Р, СА с, 


пах 2 
(4n) IN А 


Since in a pulse radar the average power is 


ВЕ. = Е 
т = pulse width 
ров] есстоп Lrequency 


and the receiver bandwidth En is approximately 1/t 


PEAGA? g 
p^ = ам Е (3-7) 
А вЕЕ (м). 
Iu o min 
ubere (S /N ).. . is the minimum signal to noise ratio at the 
о отіп 


mut Of the linear section of the receiver, necessary for 
ШОШ лоп. The first important characteristic in equation 
EX iS the fact that some of the variables can be considered 
deterministic but others must be analyzed as random variables. 
Nx ndom variables are c, and S/N. Statistical discriptors 


D 


for o, are Mendy extensively studied [1-2] not only in the 
ШЕШ Та! case but also for some specific cases, as for different 
ве ОЕ сосет. Since the noise is normally treated as 


stan, all necessary statistical discriptors are available. 


iS 





EDS рано раса бтссттрсотв together with the decision 
teria dc termine the probability of detection and 
meobability of false alarm. 
B. RELATION BETWEEN THE RADAR RANGE EQUATION 
AND SIGNAL PROCESSING 
After the War (1945) most of the effort in increasing 
ее capabilities (eq. $ 5-7) was focused on increasing 
average power Ру!» antenna gains (G), and mixers and 
receivers with lower noise figures Сех ПИР А Пе о пе слона 
che klystron, the peak powers reached the megawatts range. 
ШОО iat was Verified that since the technology was already in 
the limits of transmitted peak power, any increase in peak 
Bewer was questionable from a technical and a financial aspect. 
ESO. since 
R = x p 1/4 Дю = constant 
in R - In k * j in P 


ав _ 


ар 
R p 


1 
>| + 


That means, that the percent gain in range is only 1/4 of the 
EN cnt gain in power. The pulse width could not be increased 
much more because of the required range resolution. This led 
to the development of what today is called signal processing 
techniques. Waveforms were studied in order to increase the 
average power but maintaining rangeresolution. The basic 
faenmrques will be presented in Chapter IV. 

The brute force approach to increase the range of detec- 
tion (3-7) was used until 1955. The technology then shifted 


tO an improvement of the signal to noise ratio. If the new 
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ГЕ to П0156 ratio is 
(S/N) = м, ШІК 


ШОО О ОУ there 1S an increase in the maximu.. ranga of 


Моресбсіоп (2-7). 


ISE ОЕ СОМРОТЕКЗ АЗ SIGNAL PROCESSORS IN RADAR 

ШЕуапрс апа ахітпеі resolution became finer, as multi- 
Beer processing due to digital techniques became more 
effective, the conventional final processor in the radar 
system, the man, gradually was substituted by the computer. 
The total amount of information and the necessary speed for 
messing, tax the capabilities of man. So, digital tech- 
memes are gradually being substituted for not only those 
processors that used analog devices, but also where human 
@eeratOrs were used. 

ОС Gives an Overview of the techniques presently 
merece in digital signal processing. As in the analog theory 
вного approximation is used, So, using the theory of linear 
EM rete time invariant systems there are two common approaches: 
digital filtering and spectral analysis. Using the theory of 
Gaetal filtering, algorithms can be implemented that correspond 
(псе impulse response filters (FIR) or infinite impulse 
response filters (IIF). Using spectral analysis, two methods 
are possible: implementation of Fast Fourier Transform (FFT) 
meeeorichms Or use of statistical spectrum analysis. 

Poridve@taletiiters and spectral analysis are affected 


by quantization noise, which must be taken into account. 
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EN "the radar appears as one of the prime 
бр cations Tor digital systems. 

uu cnn ous that the first low data rate digital 
processors were implemented in radar systems, basically to 
perform tracking and weapon order computations. The next 
Step was the implementation of the MTI processor. А flexible 
E cammable digital processor was not achieved until very 
Mmeeemtiy. The low cost of semi-conductor memory and the 
advanced technology in microprogramming, with inherent increase 
ШІ өресі, were the main reasons for this achievement. Now а 
Васа] cencral purpose signal processor for radar will be 
described. Figure 2 shows a block diagram of a digital signal 
Processing system for a radar. The main processing unit is 
Ma rithmetic pipeline, which is a general purpose signal 
processor (GPSP). This consists of five sections selected 
in a way that the most frequently necessary operations can 
mempertormed; Fast Fourier Transform (FFT), recursive and 
 Ситвіус digital filtering such as MTI, threshold gener- 
memoir, реак detection noncoherent video integration and range 
Ine le estimation. The first block in the pipeline, the 
Matrix switch, selects data from data memories, also from the 
uc cr through a track buffer, or from anoth er GPSP via a 
Dx u5und connection. The second block, data scaling, 
ние (пс data for floating point arithmetic. The complex 
Nip ication block does the equivalent for four real multi- 
plications and two real additions, which corresponds to a 


we multiplication. Since some operations, like integration, 
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The re-normalization block is sometimes needed. Connected 

to the arithmetic pipeline is a memory section and a control 
Eon. the output of the arithmetic pipeline is relayed 

EN unput/outputocontrol unit (IOC) which is a buffer for 
ШӘ тайат data as well as a macro control for the processing 
lNxUgons ot the GPSP. The advantages of this processing system 
МПа: олсе the data is started the process is continuous, 
amdedoes not follow the normal sequence of fetch, operate, and 
Eu characteristic of general purpose computers; the pro- 
mor has a horizontal structure as opposed to the normal 
Па! structure of computers; this means that an instruc- 
matter execution, causes all operations specified to be 

ФР utedgd in a clock cycle; there is a hardware separation of 
NS IUctions and data in the memory, which simplifies the job 
Ae programmer; the control memory is made oí RAM's which 
ие a much greater flexibility in subroutine changes; and 
Mmemexecution of instructions is within the same time frame 
ШЕ соттесропав to a range or doppler cell. 

M O eS interact with the arithmetic pipeline. Two 
em store and recirculate data, the third stores weighting 
Ments lor such filters as MTI, FFT or pulse compression, 
OUT tn stores micro opcodes. Data sequences and arithmetic 
n rons-are both controlled through software, by two control 
ШЕ he program that resides in the sequence control memory 
(SEQ-CONTMEM) , informs the machine of range and doppler dimen- 
ЕП ОТ the problem. The program in the Arithmetic Control 
Memory, ACM, informs the machine of the particular algorithm 
meme Used. 
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ВТО ASPECTS OF THE THEORY 
ОШ АШАК ЫЛ ОМА PROCESSING 


AS TGNAL PROCESSING 

The desire to transmit information is closely related to 
the processing of signals. Since normally the information to 
Dewtransmitted cannot flow through the available channels, it 
is necessary that it be processed to use the chosen channel. 
Moon? at Fig. 5, the inverse process has to be accomplished 
Oder to get the information in an understandable form at 
 тгессіуіпс end. The process is not so simple, since in 
general the processors themselves and the channel introduce 
nosse that may or may not eliminate or change the understanding 
e messege. So there is a need to interpret the results 
some decision criteria. Here decision theory plays ап 
portant role. Even for the same type of application, for 
Mitweance, radar, decision criteria are not the same. The 
Euwerra may be of constant false alarm rate, or may be of 
INS probability of detection. Basically every application 
EE stem must have its decision criteria. 

mine are two main approaches to the processing of signals, 
the time domain approach and the frequency domain approach. 
LI rhe fact that they use completely different hardware, 
in the pure mathematical way, they areduals of each other 
Simeon the Fourier transform pair. So the only reason for 
EM (e OL approach is due to more perfect reliable type of 


Hardware. 
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Ma mera dvance ot computer techniques, both approaches 
are split between digital and analog signal processing. 

It was only after World War II that electronic engineers 
became interested in the applicability of digital hardware 
techniques in signal processing areas. But it was not before 
Aeae sixties that digital techniques suffered a big jump 
momeacing some of the analog techniques. The foundations of 
ИВ: 1 signal processing can be related to the Laplace's 
 ітеіоттп theory; but only in the mid sixties was the theory 
mu cented in a formal way, Various papers appeared at that 
related to this subject, but only in 1969 was the first 
Sot made for a comprehensive theory of digital signal 
доо по [21]. Most recently (1975), two very comprehensive 
books by Oppenheim [3] and Rabiner [4) can be considered as 
Eu can excellent treatment of the subject, not only through 
meemmethematical structure but with very good applications. 
Today's trend is definitely towards the substitution of analog 
шишево5зіло by digital processing. The speed achieved in 
Eus digital processors, the very fast algorithms used to 
ENENent the FFI, the gradual substitution of infinite impulse 
ие methods to finite impulse methods due to a better 
Ее of theory and higher efficiency of calculations are 
ШИЕ or factors responsible for the increasing shift to 
auc cechniques. There are yet some type of applications 
Du he speeds are so high, or the digital hardware is so 


сатса that analog techniques are still used. 
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B. AMBIGUITY FUNCTION AND AUTOCORRELATION 
У сли Росо 

Considering only non-accelerating targets, the two 
ds or imtormation of interest are the position of the 
Em rand ts velocity. Given two targets, the capacity 
to differentiate between their positions and velocities is 
DN for most radar applications. 

ИМА о to dırterentiate in range is called 
MENE” resolution, and in velocity, velocity resolution. 

СИРЕ e resolution is concerned, it. is 
ИЕ the shorter the pulse the better is the differentia- 
И реткеепл the two targets, that is, the higher is the 
mange resolution. With the measurement of velocity, since 
NEC dectly related to the phase difference between sig- 
РЕЙ е longer the pulse the greater the number of cycles 
Сап be compared and consequently the better the doppler 
Ie velocity resolution. So, at first glance it looks 
EN No optimize one of the resolutions the other will have 
Паскал Е1сеа. The other solution is to try to get an 
Bm solution for both cases at the same time. There is 
ПИТИ Ра сі уе мау to express these conditions in a precise 
мау; that was the reason for the appearance of the mathe- 
MC oncept of ambiguity function. 

Range Ambiguity Function 

Бес the transmitted signal be represented by S(t) = 

Bee), vlt) = u(t) exp. ) og t. Then, with no doppler 


mulus tie received signals from two different targets are 
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ПОВИ СС о) 


I ОООО me difference between the two targets. 


If it is chosen, as a measure of the resolution, the mean square 


of the difference between the two received signals, then 


Hc Ue - v(t-)|*dt 
ШЕСІ ас. HIC к 
Р (t= 1) w*oE)v(t- vat 
-F oo 


Ze) de = 2Re[exp jug s[/u* (t)u(t - 1)dt]] 


cer the first part of the result is proportional to the 
total energy of the signal and has a constant value, only the 


second term is going to make e(t) vary. Therefore 


+ 00 


х = Е 
HRS CEJUt > T)Jdt, complex 


conjugate 
ШШегіпес ас the range ambiguity function; only c(t) varies 
Кила је о: e£^(t). As can be seen, there is a perfect 
ШҮЙСІТУ between c(t) and the autocorrelation of the signal. 
The ideal case would be to have the ambiguity function with 
a Spike at the origin and zero anywhere else. That would mean 
mc only situation where it was impossible to differen- 
IEuEEthe two signals was when they were on top of each other. 
As an example. Fig. 4(a) shows the autocorrelation 
Bon ot a pulse of duration T, Fig. 4(b) its uncertainty 
ние о (т) |". In part (c) there is a rectangle with the 


same height of c?(o) and with the same area under |с(т)|“, 
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(a) 





Су) | 


(р) 








EIE 1. 0 4s ft can be seen, the base of that rectangle 
15 a measure of the spread of the curve and is called the 


delay resolution constant; so 


m - аат 


е2 (о) 


 Гітс5 1s more practical to define Ar as a function of the 
+ оо 


Каса Let u(w) = Flu(t)] = f, u(t) exp - jut. 


since by Parsevall's Theorem F [c(t)] = |u(w) |? 


гп Ј  |и(о) | "да 
AT = 


[/ 7 ]u(u) fd»]* 


тире effective bandwidth is defined as 


pedo] 





W A ss v UE a 
e ies | 
471/.. |u(o)|'do 
Ат = 1 
2 W 
e 


Ehe range resolution becomes 


CAT 
AR = Es 
these definitions make sense since as pointed out 
before, with a short pulse there is a higher ran® resolution, 
but to have a short pulse there is need of a wider bandwidth, 
НИШ 15 саилуајепсе in the two expressions for Ат. 
eee Velocity Ambiguity Function 


Using the same type of mathematical logic and defining 


ФО) = Flo(t)] 
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ГІШ ех correlation function is determined [5] to be 


H 


К(Е Ла* (21 Р)и(21Е - 2тЕРаЕ 


Ер, 


МА) | ~ ехр јопе tdt 


а 


Acne doppler resolution constant [5] 


Af, = ЕЕ ad Aut) de FL 


К? (о) ПА re 


Шо Те 15 the effective duration. 


2 = 210 
C 
ИН 24770 
dfg = Eo 
which infers Afq - ¿2 av 


mie velocity resolution constant is 


Meemeraphical interpretation for Te is similar to the 
пса interpretation given for At in Fig. 4 but instead 
EN uns the function |c(1)|*, |K(fg)|* is used. But as seen 
EN T. the optimization of one resolution implies the mini- 
О ion ер ет ле implies the need to study a two 
dimensional correlation function in order to see the mutual 
effects between the two resolution parameters. 


DO if the transmitted signal is 
w(t) = u(t) exp j2nfot 
mit-m@eecived Signal with doppler and delay will be 


Пе t-1) exp j2n(f, - fq) (t-tT) 
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The mean square of the difference 


є 2 


f|w(t) - v(t-x)|*at 


Утес а two dimensional correlation function of the form 


Ir. fa) ПОО Ет) exp-J2mfgtdt 


Er 3t can be seen that 


mtr, 0) = c(t) 


RO, fa) = КСЕ.) 
If the volume under |Х(т,Ғд) |2 15 determined and 
divided by |X(0,0)|*^, an equivalent resolution parameter 
Gaptea the effective area of ambiguity is obtained. 


aC = Аха Она — 


|Х (0,0) | 2 
[ers easy to prove [6] that 
Nu. £q)|?drdfq - |X(0,0)|? 
Ет 


So 1ї an effort is made to improve one dimension, the other 
never get better. But given one acceptable resolution 
Ne. dimension, there are ways to find the optimal solution 
meee Other dimension if X(1,fg¿) is known. Fig. 5 is an 
ambiguity surface for a monochromatic pulse with а smooth 
EE. It would be a reasonable temptation to try to 
design the waveform from a desired ambiguity function. That 
BEN out to be impractical. First, a criteria to design the 


BE tv function would be difficult to build. Secondly, and 
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most importantly, when designing a radar system normally there 
Ex other more important factors, from economic to space and 
ШОЛЕ constraints as well as technical, that dictate the 
boundaries for the waveform that we have to pick. Also, there 
15 an important relation between the chosen waveform and the 
mae Ot clutter model, and the approximation of the clutter 
model. Normally some flexibility in the digital processing 


is desirable. 


EN PULSE COMPRESSION 
I Concept 

Soon it was found that trying to increase the maximum 
range by increasing power, keeping the same pulse width, was 
meesdafficult way. Isolation problems were difficult to 
handle as well as reliable components. The only way to in- 
e the energy of the signal was by increasing the pulse 
ии. But an increase in pulse width would reduce the range 
ШОО топ. Since the real problem was to increase the time 
of transmission but at the same time do not decrease the 
bandwidth (effective) of the waveform some different types of 
waveforms were studied, using the ambiguity function concept. 
NU ocessing of the signal should be such that the trans- 
mitted waveform with high time of transmission and small 
Ee Л уе bandwidth (We) should be converted to a high energy 
e ith pulse width approximately 1/W,. This is called 
pulse compression. There are two basic ways of implementing 
Se compression techniques. The active way, which is a 


time domain approach which basically uses correlators in the 
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@eecction and active elements to produce the waveform in the 
transmission. The passive way, which is a frequency domain 
approach, uses passive filters to generate the waveform and 
matched filters in the receiver. Each of the two has some 
EERUStrons, and it is even possible to build a system as a 
combination of the two processes. Table I gives a relative 
performance of various types of pulse compression techniques. 
E Linear FM (Chirp) Pulses Compression 

er mNeMPUSSSTVe generation ef Linear FM signals 
fier be analyzed. An IF pulse generator feeds a dispersive 
Гау filter with a frequency versus time characteristic as 
EU v. O(c). The signal is then up converted and only one 
(Пе bands is transmitted. At the receiver the signal is 
down converted to JF and then passed through a dispersive 
Шин with the opposite slope (Fig. 6€ c). Using at recep- 
tron a mixer sideband inverter, the same type of dispersive 
filter can be used because as it can be seen, Ё, апа fz 
о Б) have opposite frequency versus time character- 
Bear, In fact, with this technique the same dispersive 
Ме for both transmission and reception can be used. 

Me analysis approach to this type of processing Is 
ео Том пс: 


Let the received signal be of the form 


КОЕ огоо сое e 2ut*), |t| « T/2 
NUS t T/2 
where u = il A eoun Slope of the dispersive filter. 
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T = duration of the transmit pulse envelope 


за = doppler shift. 


The pulse spectrum is 


+ со 


F(w) = E AA 
+ 00 1 
=). coUo t са O t zut ] dt 
walter transfer function is 


Шо) = ехр ј[ (аб - w)^/2u] 


which implies G(w) = H(w)F(w) = exp j[(wo - w)?/2y] 


12 1 
j | exp j[(wo + wg - o)t * 5ut*] dt 
2 


ИП ітріісс the time function at the output of the filter 
15 
_ 7-1 ie NR 
e(t) = F " [G(u)] = т 7.60») Xp ots. 
+ оо 


which implies g(t) - d 4, exp j(ug - ш) 2/20 


ue 


Т 
Ке РА Г а T а 
-T/2 


Е ЗЕ Ча 


Imverting the order of integration: 


T? + 00 2 
| 1 ОСО, 
Í С ЕЕ E + Core) сс 


Es 


dwdt 


g(t) 


IA 1 
ето хо іо т оат + vs + 


so | 
E72 H 


X 


+ оо ТЕ ш | 
S Exp 1 le Mess 25 m = P алш ат 
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поинт = nt) 
ма 


E Wl v 


2 
же UO x n 2 + Е І 2: _ 
V 2y 7 MT 007" > ШЕ UTE Wot 


Multiplying g(t) by exp -v? and exp v? 


I | Т 
met) - Ау хр j [vot % шат % > Wee = РУ 


X J exp j [w - Ұ2ы у|2/2ң «шат 


= 2 2 2 
Since, с /2н ҮІ? у . 22 vu 
ZU ZU Zu 
2 
1 
- 2m ES u | ete) 
2 
= y? -wt - 29 y 
If u - 5 ny | | | 
E --- , (һе second integral will be of the form 
VEU 


Y2yu f exp dedu mm ехр 11/4 


li 


meh implies g(t) ars CoU) itu 7 ut? + ДӘ 


(о) 

тт aa Ned. 
Бе шот + war + > По в ати E ut? + wot 
which implies g(t) = (ay Бе = Е”) 6 ways? 
which implies g(t) = "ШТ соз (ОСЕК: > Шел OS 


Х 


where Е wt) lf 2 
e шаг Au 
X 2 + 5 t 


ББ 





Sain X 


Ben these results it can be seen that if за 7 0, а x 


meee ОС envelope is obtained with a peak at the origin of 


magnitude АТ) moo the higher AuT the better the output 
" 
ПИК ft ша 7 0 there is a shift in the sin x/Xx curve which 


implies a range error for doppler frequencies greater than 
Zero. 
meee Matched-Filter Approach 
The matched-filter approach is another passive method. 
ШИЕ analysis of the former will give us a basis for comparison. 
Ber the received signal be 


fet) = cos [(u¿t+ шд)т + 5 ne |, |6 е тулке ссн 


Me output of the matched filter will be of the form 


g(t) = Ba Nar - t) dt 


where h(t) =Kf(-t) and F[h(t)] = H(w) is the filter transfer 
топ. K is a constant factor to give a unity gain. Іп 


this case, К = у (20, 50 


Ti 


со 1 
в(тьма)з уд fo соз (шо кад б т 5 ме"! 


cos [ug(t - t) - и (т - T)?]dt 


after some steps [7] 


TUS) eV-Zw.cos[(w, + 261] sin EL (T - |«|)] 
g(t,ug ea X 2 7 





за + ут 


onsec «T 
Вт, ша) = E 20) соз [ (шо + 50) т] => , Ice 
Where х = е aye т) 
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CES NOt a linear function of т, the shape of 
шиоюсогує, unlike ın the previous case, is not a sin x/x type 
EN vc. Also, since wat does not represent a linear shift 
MENE axis due to the non-linearity of v with respect to x, 
но 7 0 the curves become distorted and not symetric with 
NI ect to the vertical axis as in the previous case (Fig. 7). 
Ime advantage is of course the maximization of the signal to 
Moise ratio. 

Ie mot Dsscrete Erequency Sequences in Pulse Compression 

Меше лоша sequence of frequencies each with a pulse 
Eon ОЕ т іп order to get a good autocorrelation function 
15 а very common method. 

I mdousmeombanations. from linear stepped frequency 
ME maomly chosen frequencies, are used in order to increase 
пап lobe with respect to the side lobes in the ambiguity 


fection. The basic format is that of Fig. 8 





р о 


Where the waveform of the nth segment can be expressed as 


v, Ct) = A exp j 2n(f t + ġa) 


ЭМ 









ооцЈипј 
цоцездігоз0упе 
дојолиз 


| м 
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WEP 
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CS | 
NS 


01 =/V.J 
оочу 
ооцеја јозојпе 


JIi -PINEN 


vo = y/o 


OT = My/Pm 
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ВР TE 


The mathematical analysis for the general case is 
Eu c urfricult and of little interest, Since only a few 
ШЕ are Of interest, consider the linear stepped frequency, 


meee is, a time function of the form, 


De 
v(t) = m na) - u(t-(n+ 1)1)] cos Су ELE 
where Dos lowest frequency to be transmitted 


Aw = frequency spacing 


N = number of frequencies in transmission. 


повен Chere 1S coherency between all frequen- 
mes, that is, o, = (ЛР оте пасош а о реге рате using a зва шоке 
с брезаљег with a master oscillator, the matched filter would 


memes) the form indicated in Fig. 9. 


ms Eb - -- ЈЕ 


Ж 
A] А 


СОЕ СЕСЕ S [ОГ 


output 
Ја то 
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I o ctherc will be significant range ambigui- 
BL ALSO, rf tAf < 1 the subpulse filters f of Fig. 9 
foe Overlap. That is the reason why in most cases TAf is 
EE qual to unity, that is, т = iT With that assumption, 
mer convolution of the impulse response of the matched filter 


ШЕ! - v(-t) will yield an output 
¡NA 


E) - Еи е ехр ј (зо + пло) Е 


мћете ШЕЛІ < мт. 


t А, ә сот all n, then 
М-1 


eg t) - exp j wt nèg (exp DE 
Since lexp j Aot] « 1, then 
eg (t) - [exp ПОБИО ы С) x sin N (Aw/2)t 
2 ZE 


ПРВО О атол the autocorrelation of the rectangular 
EE. the final output time function is 


Et) s (1- [t]Af) sin N(Aw/2)t 


(N-1) Aw 
SN 2 


А | 
COS А lt 


RO shows the effects of tAf # 1. The nules of 
the output waveform occur at sin NAw/2t = 0. 


Which implies N вы Е = + мт 


не 50, the first side lobe is 13.46db below the main 


lobe. 
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SUBPULSE AUTOCORRELATION 
ТАГУ ìl 








Percy important quantitative factor in all these 
waveforms with sharp autocorrelation function is the pulse 
Bempression ratio. It is defined as the ratio of peak power 
after compression to peak power before compression. Table II 
Res some data including compression ratios for some typical 
EU cr pre-pulse generation and compression circuits. 

Another important quantitative factor is the side lobe 
meauction factor which is defined as the power ratio of the 
fem Lobe with reference to any side lobe. 

Pependine con the specific application, it is sometimes 
preferable to sacrifice the main lobe peak power or even width 
КОО Б point), but get a higher side lobe reduction. In the 
radar problem those side lobes may be identified with targets, 
the need to suppress them if wanting to use a larger dynamic 
ШЕШЕ 15 obvious. 

There are various types of weighting functions. Тһе 
weleghting is applied in the amplitude response of the matched 
вет Jn Order to deliberately produce a mismatch. 

ШОШ requeney response of the matched filter is such 
ШИЕ (Пе signal to noise ratio is maximized. So any mismatch 
will lower that signal to noise ratio. 

Mito ss in signal to noise ratio duc to weighting 


Шо еса the loss factor І It is defined as 


2 
[S (S/N) weighted _ [/ то (Е) dt] 
р (S/N) matched EU uode 


w(t) 


powchtang function 


Т = processing time interval or the length of the 


eane mitted waveform. 
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lU case of discrete frequency amplitude weighting, the 


INT. factor is 


N 
ee 
L > > [8] 
N [z А21 


ШЕ ртор1еп оі меісһібіпр 15 in effect a way of producing a 
pulse compressed spectrum that yields the wanted waveform. 

In the continuous linear FM case, the spectrum of the 
Euch filter is of course a pulse in the frequency domain. 
EE т һе inverse transform of a pulse is a sin x/x type of 
Bometıon, that is why for the linear FM waveform a sin x/x 
Best output is obtained. From the Fourier transform pair 
it is known that a wide pulse in one domain corresponds to а 
EU pulse in the opposite domain. So the amplitude weight- 
EN UuD»ctions that will lower the side lobes will have a bell 
E Dered shape characteristic in the frequency domain. That 
EN the time output of a matched filter with a pulse in the 


frequency domain as transfer function, 15 


Bl. м : 
et) зо Л, a(w) exp j wt dw 


where a(w) is the weighting function; the most commonly used 
рон по functions are: 


NL Ncos)ine function where 


a(w) = cos p терен а (оо Isat the cemer of не 
рес ктш 
© (£) - ТА a п 2 COS Ut 
©, eg 25 Le 5 жу exp J w E 
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The main lobe will be between nules. The nules will 


EScur at cos Wt = 0 


which implies wt Е т ТЛ ТЕЛИ 0% 
PN 


Or О ОВЕ 
ZW 








Considering the lobes between nules the side lobe reduction 
5. 5 db. 
Another class of tapers are the Hamming functions 


which have the general form 
СОГ о = a + (1 =a) cos (2) qs T 


Es the center of the spectrum is assumed to be at о = 0. 


In this case the output will be 


E = mac t | ze 
po P6 ШІ Zu It) 
КИО = .54 the side lobe reduction is 42.8db which represents 


a much larger dynamic range. Table III gives data on side 


lobe suppression on frequency coded waveforns. 


D. MOVING TARGET INDICATOR (MTI) 






АЕ а high pass 









| of a variety of types of clutter which mak: 


mer, If only the elimination of one type of 
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effi ciency aot 
the same processor, the efficiency drops. 

Pepending on the parameter uscd, the MTI processors are 
Шивисаїіу divided into three types: those that use the phase 
information of the returned signal, those that use the ampli- 
me, and those that use both phase and amplitude. The basic 
mera curation of the phase processing or coherent MTI is that 
Мао. 11. The use of two oscillators, the stalo and coho, 
Is for up and down converting the waveforms in order to get 
BE (rect phase comparison, since at low frequencies the phase 
асу measurement increases.  Coherency is also obtained. 

Consider two signal returns from the same moving target 
ІП а separation in time by T = 1/fy o рые 


R 
НС топ frequency. The signals presented to the summer are 


Ber E sin (шас + $4) 


t] 
H 


E sin (ма (t * T) * фа 
41 7 b> = E 
lu Noutput of the canceller wili be 


| ЕБ 
И 2 E sin ES Gos Ga] | 


Det is seen that the cosine waveform is modulated by а 


sin ES . So, independently of t, there arc some doppler 


frequencies where the power out will be zero. These are e 
med luencies where Е nie „Пар то» Ра = nf, Е Ве 
ÍR 


(ПС рпісс repetition frequency). 
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EQUALIZ, 


E, = E(t)-E (e + 1) 





E (1 + Т) 


48 





pace the doppler frequency is a function of the radial 


2V 
a locity Ур; - = LA meciere are also blind=specds tat 


correspond to the blind frequencies (Ра = лї) 


Af 
V lind = ul oS s ап antecer: 


Pie. 12 shows the power response for a single delay line 
р тосе$50от. As it is seen, the multiples of the PRE are 
Ки е5 but the other frequencies are not processed in the 
Same way. That is very far away from the ideal case. 

mre Clutter spectrum beating with the pulse repetition 
ШОЛ епсу 15 translated to every multiple of the PRF. So if 
Mmmemelutter spectrum is, after translated, of the form indi- 
ратес іп Fig. 15 (a) the ideal MTI processor should have a 
Mmeeponse as in Fig. 13 (b). 

МИТЕ 15 known the single delay line processor response 
a away from that of Fig. 13 (b). One way to improve the 
EM  ocessor is by using multiple canceller filters. 

lan n-stapge cascaded canceller (Fig. 14) 


E 
69, - sin ES [9] 


where SR and 5. are the peak output power and peak input power. 


ле ртођјел 15 that if many cascaded blocks are used, the 


2n 
(sin 119 Auction may cut some low velocities or interest. 
R 


ie transmitted waveform has a spectrum as in Fig. 1. 
The fact that there is more than one line in the spectrum 1s 
Bus importance to the MTI processor since all the lines are 


at blind frequencies. 
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ШИЕ ашрілсійе processor, also called noncoherent processor, 
has a block diagram as in Fig. 16. The advantage of this 
processor is due to the fact that stability problems almost 
EM uppear since there is no coherency. Since most high power 
Mmeansmitters are not very stable, this is a real advantage. 
ШЕР ріс disadvantage is that it needs a permanent existence 
EN strong clutter return since without it, it won't work. 
17 represents the phasor diagram that explains how the 
memeoherent processor works. Since the clutter is of the 
Form, p ee zn ІП а рлазог атасгаш 15 represented бор а 
СЕЗ олату vector. 


Bet E. (t) be the return of a moving target. 


Then E. 

E, (t) = | ES Ct) | cos wat 

E. (t +T) = Е. СЕ) | COS (ша С O] 
Bet Wat = > - Аф/ 2 


! 


ша d) ф + Аф/2 = ф - 9/2 + wal 


Domeh implies ша! Аф 


Ме кота] signal at time t (Fig. 17 b) be 


-> -> -> 
ES E. + Et) 


ПОТЕ total signal at time (t * T) (Fig. 17 c) be 


> yd > > 
EEUU E(t + ВЕ Се jm 


From Fig. 17 (d) it can be seen that due to p there exists 
An amplitude difference between [E] (t) | and [Е (Е). This 


difference in amplitude is going to be used to identify the 
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(а) CEM) = E(t +7) 








(b) 
(С) 
Е,Е-ІЕ| 
2 
poc» 
Е. 
(а) 





moving target. In order to get better results, instead of 
the difference in amplitude, the difference of the square of 
the amplitudes will be used. That is the reason for a square 
Шек detector in Fig. 16. 


Rom Fig. 17 (b): 
E = Be + В. 1 26782 со$ (ф - дф/2) 


(the equation relates only the absolute values of 


the vectors) 


шоп Fig. 17 (с): 


ПБ 2E.E, cos (d + A6/2) 


2 
2 
(only amplitude relations) 


Then E? = E? - E2 - 4E E. sin $ sin 2£ 
1 С © 2 


Í 2 
since ф = за? + Ag/2 
and Аф = шат = 2744 т = 2 
К К 


which implies 


2 = . T f . 
Er AE.E, sin F Cae (eet ey 2) 


50, the blind speeds are the same as in the coherent 
processor, and the transfer function must have the same shape. 
But if Ес = О that implies Ec = 0. That is the reason why 
Eu lytude processor must have a clutter return, and the 
stronger the clutter return the better since Е. increases 
linearly with Б. 

mkt ter cancellation can also be made at IF frequencies 


ВОЕН of at the video part of the receiver (Fig. 18). 
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Ee two lF signals separated by Т = ES be 


ER 
E, З E sin [27 (fig == tt > Фо] 
Ба = E sin [2n (fi, + f4) (t + T) - $5] 
where Фо = о 
C 


Б. = Ei - E, ТИ т (р ЁТ] GOs 


[2n(fgp* £D (t* 29 - Фф] 


che output of the summer. 


The output of the phase detector would be 


ieee ee san [nck 


5 EQ)T] eos i2m[fjt* (fry f) 2] "M 


+ 
IF- О 


Еше the video envelope is of the form E sin [1 (fip + Êa) T]. 
Се іп ап MTI processor if fa - 0 the output must be zero 
that means that ПЕТ ас ЗЕБ Лат Ето те. 
fip = T n fy, pcd freguency must beaa multiple- of the 
Meese repetition frequency. 

mmemextistence of blind speeds can constitute a problem 
Кише еџету velocity is not processed in the same way. One 
Eee most often used methods to improve that situai rongis 
ШОО е use of staggered PRF systems. There are some quanti- 
Пи гастогт5 that describe the efficiency of an MTI processor. 
 ПСС (Ппеу аге the descriptors of the MTI processor, it is 
Паје to state their proper definitions. 

Improvement factor is defined as I - io where rj is 


Miewoucput target to clutter ratio and г; the input target 
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| 


ИО СТ ratio. This definition reflects the gain as well 
EN Uc clutter rejection of the processor. 

шесте Visibility is defined as the capacity of a 
ша» to detect moving targets in a clutter environment. А 
Eur with x db of SCV is one that is able to detect a target 
Ea clutter that has a signal x db stronger than the moving 
@ereget. The SCV cannot be used as a parameter for comparison 
Een radars since the target to clutter ratio is a function 
Шоле size of the radar resolution cell. 

ШЕЕ two main problems with MTI processors design are: 
SE, a correct modeling of the expected type of clutter, 
As econd, the correct shaping of the transfer function 
аеш по оп the type of clutter. Іп [2] there is an extensive 
@eecription of the different types of clutter and modeling 
processes. | 

Eun the platform of the radar is moving, the stationary 
targets will have a doppler return different from zero. This 
EN typical case of the Airborne MTI processors. Various 
 Пүйиес оғ clutter looking and automatic tracking must be 
Пи ihe techniques get much more complicated and the 


processor becomes much more dependent for specific application. 


ПОС ТЛІ SIGNAL PROCESSING, DFT AND FFT 

Due to high speed, low cost, versatility and almost inde- 
КООШО С Of external condition of recent digital computers, 
ENDE eessing of radar signals in digital form covers now a 


wide spectrum of the radar signal processing. 
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The two basic methods of approach are represented in 
En x diagram in Fig. 19. In the first case (Fig. 19 a) the 
Signal is first sampled and quantized in such a way that the 
Signal is converted into a sequence of numbers. The numbers 
ле леп stored and the processor, using arithmatic and logical 
Mmerations, manipulates the numbers according to an algorithm 
which is a function of the type of filter we want. The 
resultant numbers are then dequantized and passed through a 
Вора filter (D/A). The second method (Fig. 19 b) is 
basically a frequency domain approach of the problem. The 
Eun tS quantized and a discrete Fourier transform algorithm 
EN  Dred to the quantized signal, The spectrum of the signal 
тел weicshted as a function of the type of filtering re- 
нирок» Then, the inverse operations take place. Due to 
NH aämprovements in faster algorithms to implement the БЕТ 
this second method is becoming widely used. The theory behind 
NIMUSt method is based on the Z transform theory [10]. The 
EUM ble Z represents a delay and is defined through its Laplace 
transform by 


Жа ехр 51 


where 1/T is the sampling frequency. Given an analog function 
NESEUthe transform of f(t) will be 
== $ n 
СА) 120 БТР 
ПИО basic configurations are recursive and nonrecursive 


meer. In the nonrecursive case the output is of the form 


й | 1 
ПУА 


-n 
O 1 ЛЕ 


60 












540)06б6! HUIUDIOM 
MEM USD > 





JO 
idtm 





(0) 





indui 
бојеџе 





ino | Јо! 


Борио пед бІр 





61 


e a танро но e EI PI eR ЮТЬ о оо рань = ние ми ар MES eR ye аа ана о жк а t, n dm amm s ap Ham s rr ee a E па — - 





штете the coefficients А. arc Obtained Irom ihe coeli iee 


of a Fourier series expansion of the frequency domain function 
Ehe desired filter. 

The feedback or recursive filters are normally obtained 
Eus the analog transfer function of the required filter and 
puUestituting S by a chosen relation between S and Z. Usually 
the bilinear transformation 


wen 


1 + 71 


15 used since with the correct sampling frequency it is very 


[11] S = 


| dm 


accurate and easy to implement. Figure 20 shows the basic 
configurations for the nonrecursive (Fig. 20 a) and recursive 
ШИ 20 b) filters. 

Eon radar applications the discrete Fourier transform 
EN лр used especially as a bank of filters for obtaining 
Aer outputs. The understanding of what it is and its 
limitations are thus fundamental. 

Ші отдет to get a Fourier transform with a digital pro- 
Gessor not only the input waveform must be in digital form 
but also the values of the Fourier transform must be repre- 
sented by discrete values. Figure 21 gives us the graphical 
Шта топ об the discrete Fourier transform pair. The input 
Burm h(t) (Fig. 21 a) must be sampled in order to be 
manipulated by a digital processor. Тһе sampling operation 
ШОО same as multiplying h(t) Бу а string of delta functions 
EN 1127 b). Since multiplication in time is the same as 
convolution in the frequency domain, the resulting spectrum 


EN t of Fig. 21 с). Here it is secn that if the samples 
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ame not at a frequency higher than the Nyquist rate there 
will be an aliasing effect that will distort the waveform. 
ME since in most systems, and especially in any real time 
radar system, the waveform must be processed in a finite time, 
Bere are for computation only a finite set of samples. This 
equivalent to multiplying h(t)A,(t) Пе A еј 
Meerangular function (Fig. 21d). Again, the previous spec- 
trum must be convolved with the spectrum of the rectangle to 
EX ue result in Fig. 21 (e). As expected, the larger the 
То? eme maller is the distortion in the spectrum. But as 
EU before, the output must be in a digital format. So it 
EN cCessary to multiply in the frequency domain [H(f) * А (1) 
Gee, by a string of delta functions, separated by 1/T,. 
Er equivalent to the convolution in time of a string of 
impulses separated by 1 with h(t)A lt) ӨСЕ Тһе лал тос 
Ето. 17 (g). Ас сап be seen, the output of the DFT 
processor is not a sampled replica of |H(f)|. Aliasing effects 
Кила се time of processing are responsible for the amount 
uU сок Стол introduced. In order to get the mathematical 
Exession for the DFT, it is only necessary to make a parallel 
Шема оп to that of the graphical derivation already dis- 
cussed. The time domain expression for the output (Fig. 17 g) 
na Obtained; then the Fourier coefficients of that 
ОШООО function represent the DFT. 

EN et) be the input function. The operation of sampling 


MS equivalent to a multiplication by 
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A (t) it á(t - KT) 


K 


б ОКК кту Б (с KT) 
К= - о 


T 
1 > Е. < Пи ерул 
0 otherwise 


SO h(t)AS(t) 


лисе ӘСІ) 


16 to be multiplied by h(t)A, (t) 


N-1 
h(t)A, (t)x (t) - и 8 (t - KT) 


where NT = de 


then we have to sample the Fourier transform of h(tJA S (t)x(t) 


Earring of impulses in the frequency domain 
Е үт ут) 
і 2 
O infers A,(t} =F [4,(f)] = T, MR Gere то) 


must be convolved with h(t)ag(t)x(t) to obtain the time domain 


answer., 
hg(t) [h(t)a, (t)x(t)] * a, (t) = 
N-1 Коо 
IE AR Is е кетт) 
К-0 
ко N-1 
= Т. о 


Since hy (t) ioe periodic function with a period NT, in 
(СО аус1і0 time aliasing x(t) is chosen in such a way that 


NN nd points of x(t) don't coincide with sampling points. 
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If they did coincide, it would generate an additive effect 
at the boundary. But NT is still equal to Ln АТО Ее 
mourier transform of a periodic function is given by the 


Mier cocfficients of a series expansion. So 





= d = s = 
Hy (nf,) Hp Gr DE 6 (f nf.) 
Т -Т/2 
where Cn = = p hy (t) exp -j zu Ute pes OSEE EI 
om 172 O 





i! -Т/2 5 ix 27nt 
EC = то / Ta оо Ро ПЕ ео) e ја: 


Since the integration is only over one period 12 





су? ГА) r 
= y ОЕ - KT) exp jet СЕ 
ET 2 К= 0 O 
N-1 еу? 
В зако ° охр- ЈЕ Е 5(t - KT) dt 
K=0 -T/2 O 
N-1 є 
= È h(KT) exp -j EA 
K=0 O 
исе То = МТ 
МЕТ 
со зо ХО В(КТ) exp -j En Ыш. пз А 


ШО к-о 
and the Fourier transform of hp (t) is 


+00 N-1 r 
E _ 2 а  Шак f odis 
Hp Ст De pos exp ЫШ NT) 
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ВЕР Be oe 
lo see that Н (дур Е песен tO 


show that C w rodic with a period NT. 


ШЕЕ п = т + М 





| | = М) _ . 21K 
Since exp -j TAE ер MA 
which implies C (r+N) = ca (r) 
: : ; т + М А ї 
which implies Hy ( NT ) = Ну (m) 


EUcre are only N distinct values that can be evaluated. 
Normally the function 
N-1 


H E = do h (KT) exp -j 2unK/N, n = 0, 1... N-1 


| 


EN ed the discrete Fourier transform (DFT) which relates 


N samples in the time domain to N samples in the frequency 


Semen. The inverse discrete Fourier transform, is, as expected 
ШЕ 
Е ml n о ie Ж. й 
h (КТ) = N c xp? ср DES О П.М 


EN uS seen that with very simple complex multiplications 
and additions there is conversion from one domain to the other. 
E east Fourier transform is an algorithm to implement 

MeS Crete Fourier transform which, by reducing the number 
нн р! зсабсіоп5 гедисе5 the processor time to perform the 
EN thm. For a better understanding of the FFT, it must be 
Шишо to a matrix form. 
Bomsiıder the DFT 


N-1 q 
x(n) = E Хо (К) exp -j a он ы ЕТ 
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ШИ! - ехр -1 Zu/N, then іп а compact form 


к= won Xo (K) 


where x(n) = (О) ‚› х(К) = хо(о) 
x(N-1) ОЕ 
yes = пи we 
е y (N-1) (N-1) 


the examination of the matrix equation shows that there will 
be № possible complex multiplications and N(N-1) complex 
additions. ТЕМ = Dale and the FFT algorithm [12] is used,the 
number of multiplications will be reduced to Ny/2 and the 
Миш of additions to Ny. Figure 22 gives a comparison 
between the FFT algorithm and the direct computation with 
NI ct to the number of required multiplications. 

Et was shown, complex notation is greatly used in the 
ШРпагіса! structure of signal processing. So the need to 
mepmesent Signals in complex notation becomes a necessity even 
in their implementation. It will be seen how a sine wave type 
EwENunDnction is converted in order to be represented by a 
Complex number. 


то. 25, let the input signal be 
EN) — A sin (wir i wit t Фо! 


where отр апа ша сап be thought as the intermediate frequency 


ENNdoDpler frequency in a radar receiver. If the local 
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Ep cullator (1/Q) generates 2 cos отр and -2 sin Шері which 


can be obtained with a 90° phase shifter, the in phase channel 


EsESnd the out of phase channel (Q) outputs are of the form 


I A cos (баб б Фо? 


0 A sin (w ¿t t Фо? 


Since the low pass filter retains only the difference frequency. 
ШОШО апа Q signals can then be thought as the real and imagin- 


part of the complex signal 
Z= 1+ j Q 
If the I and the Q channels are sampled and converted to 


sequences of numbers x, and Ук» the resultant sampled pair 


K 


can be considered a complex digital word. 





Ак E ax + j Ук > [Ах ехр ) ФК 


where x IN eels (ша КТ i Фо) 


Yn = A cos СУ $9) 


| Ак | D Фе = wg KT TED. 


assuming a sampling frequency f - n Since ФК varies linearly 
ИШЕ the complex digitized signal can be interjected as a 
Iur of amplitude A that shifts at every sampling time 


KT, by забої, 
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feo NIHETIC APERTURE RADAR 

Using the pulse compression techniques already mentioned, 
[ems possible to achieve very small values of range resolu- 
tion. In many radar applications such as aerial photography 
MemeOnly a very small range resolution is necessary but also 
ameecquivalent small aximuth resulution is desired. The con- 
ventional method of reducing azimuth resolution by very small 
beamwidths will always require big antennas which, cspecially 
wman borne radars, is a problem since space is a limiting 
Ae or. Also, the fact that azimuth resulution is an increas- 
ing function with distance makes it many orders of magnitude 
тег than the range resolution even for small distances. 

Consider a radar without a pulse compression processor. 


Then the range resolution ôR is 





ӨК = > B C - speed of propagation 
й B, 7 bandwidth 
Eee antenna has an aperture D, the beamwidth is Or = >: 
ie azimuth resolution б д7 NS 
5 EET = RA R- distance to tne target 
AZ B D 


И 


ПС B EIU MHZ, f 5000 MHZ, D = 10 yds and R = 100 mi, 


me SR = 16.4 yds 


S 17 = 1300 yds 
Synthetic aperture radar is a way of processing radar 


signals in order to get an equivalent antenna aperture much 


bigger than the real antenna to reduce 9A, . In the synthetic 


[5 





К Тиге сазе а single antenna is translated along a line. 
The received шта = are then stored, and processed after 
the radiating element Was travelled 2 fixed distance. The 
syechetic aperture Concept can be viewed through two different 
aspects. From a doppler viewpoint or from a linear array view- 
point, the basic Principle of the doppler effect viewpoint is 
that there exists a One to one correspondence between the 
long track coordinates оо reflecting Object and the instan- 
taneous doppler shift. 

Consider Fig. 24 in which an airplane flies at an altitude 
n With a constant Speed v and with a radar antenna that 
Ш патгеѕ the area A. The radar antenna does not rotate 
but simply moves With the airplane, 

If there is no pulse compression and the width of the 


transmitted pulse is T, it is known that the Slant resolution 


Constant is 
C С 
5 T 


Р: 2 Ва 





SO the corresponding ground resolution Constant is 
ор = бр sec). + > Sec y 
Also, from Pee 24 it is seen BID EE ae а 8 is the beam- 


Width of the transmitting antenna and Der is the the along 


track resolution E 


Since PES 
= = AL 
then D Log 7 Е 
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Бона о 


ИЗ 





and as pointed out, it is impossible to increase D beyond 

certain limits in order to decrease Рх- БИЕ cterclatagoH 

is found between doppler frequency and the росто Ко ЕЕ 

radar, then it is possible to make an azimuth discrimination. 
From Fig. 25 it is seen that the radial VeElOCity with 

E pect to the target is 


We ал б 
r 


Тһе doppler shift will be Ра = т. sin 0 where e is the 


@eansmitted frequency. But if 6 is small, which is the normal 


physical situation since the total illuminated area is small 


ШООК 1S big, then sin 0 = ө = Хх 00 which makes 
27 
т Ра но 


ENDS seen that with a frequency analysis of the received 
medi, it is possible to get an azimuth resolution that is 
ES function of discriminating between the closeness of 
EN requencies. The basic principle for the explanation of 
Exetic aperture radar using the linear array theory is the 
following: 

ШИ а linear array various elements are fed at the same 
time and the received echo is also received at the same time, 
Te to path differences a receiving pattern is generated. 
Meme synthetic aperture case there is only one transmitting 
Bent that moves with constant speed along a line. So a 
long antenna will be formed not by physical means but by signal 
Processing. After the radiating element has travelled a dis- 


ват се Lose Seorresponding to the illumination time of the real 


Ф 
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Шет (ле stored signals after weighting and phase shifting 
Besemple the signals received by a big linear array. 

iworGases arc important to differentiate: the unfocused 
memenetic aperture and the focused synthetic aperture. 

вене Тосизе synthetic aperture radar if ө = 0 (Fig. 26а), 
EN sena! processing (for R 7 œ) or by real situation (R = о) 
Bueesrgnals all arrive with the same phase. In a real linear 
meray the angle selectivity is provided only during the recep- 
wno the signals. But in the synthetic case, since only 
one element is radiating and is moving, it is necessary to 
E ount for the phase shift due to the URSUS d and Тесс 
ПАО, due to the fact that all the round-trip phase counts 
ШЕП (Пс formation of the receiving pattern, as opposed to the 
EMEN se in which only the received path is responsible for 
NI cceiving pattern, the effective beamwidth is 


À 
eff 


А 
ore 





instead of 


~ 


e 





| D eff L 


Все тһе length of the synthetic aperture radar is equal 


meme distance corresponding to the time of illumination, 


L = RA where D is the horizontal aperture of the physical 


eff m > 
Enna and R the distance to the target (Fig. 24), then the 


| th resolution is 


t| (C 


| The azimuth resolution is independent of A and R in the 
focused case which is of course an important result. This 


means that the smaller the antenna the higher the azimuth 


resolution. 
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In the unfocused case it is important to account for the 
phase difference between the center and the ends of the equiv- 
ment linear array. А quantitative interpretation of the 
purocused synthetic aperture radar can be explained as follows: 

Bonsıider Fig. 26 (a): Due to the fact that Ко 7 R, there 
15 a phase shift between the center and the end. The total 
NNNM sc 15 the vector sum of small discrete voltages that are 
received along the path Loff (Fig. ЛОВ) IF the Пиш оГ 
Гессе voltages is large, it is possible to make the approx- 
imation of Fig. 26 (c) and determine the relationship between 
the resultant voltage V and the linear sum of the discrete 
voltages M s MINE dD point is defined as the boundary 
шесте атор їп voltage, a must be approximately m/2. But, 
when a = 1/2, the radius ry (Fig. 26 c) is perpendicular to 


ebe tangent t implying that the phase difference between 


|2 
the voltage at the end and at the center must approximate 45°. 
Пише сет а limit on Leff’ Os ISR SO шул 


NES 71) the conditions are met. Using Fig. 26 (a): 


ВА , A? 
2 2 - 2 = 2 Er = 
ee 2)? + R (А + 1/8) R2 + = + е, 


2 
if ez is very small compared to the other terms Lee = У RA 





висе Ozap ^ „455 p^ = > "AP 
eff R 
the azimuth resolution would be Jv = КӨ = 5 УХ 
: і О ПУТ . 
Prince Önz Bu tumetion of R this Ts a Solution Detween 
the focused case and the conventional case. Figure 27 gives 


iia tion between the three cases. 
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The signal processing theory for synthetic aperture radar 
can be understood through the analysis of the ambiguity func- 
Ио: the received signal. If it is possible to separate 
вое range and azimuth as the two ambiguity factors from the 
EEnewal ambiguity function, then the determination of the 
azimuth ambiguity function and consequently the azimuth reso- 
Eun function aro immediate. 

Memsider a transmitted signal f(t). After reflection in 
HUE round or target with reflectivity p(x,y,z), the received 


бише! will be 


S(t) = 555 p(x,y,2) Е (6 - 2) ахауае [13] 


IM he integration is made over the illuminated region, 

ть (һе distance between (x,y,z) and the radar antenna. 
mremmrecessing of the signals in fact recovers a signal pro- 
ДИ to p(x,y,z). Using matched filter techniques, the 


Eee will be of the form, 


f £*(t - ARS S(t) dt 


EC 


S555 oQuys2) £ (t - 2) £* (c- 5—) dtdxdydz 


ҮШ ШЕ 15 the distance from the antenna to the point 
EN 5) corresponding to p(x',y',z'). Let's find the 


EN ey function of f(t). By definition it will be 


2 Р ZI 
ф(хру,тух' у! ,2') = f f(6 - 20) # (6 + 505) де 


ns means complex conjugate. 
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— 


— Н ЬЕ . 


50 
AI 77700 (х,у,2) dxdydz, 


is a weighted average of p(x,y,z) where the weighting function 


EN у,2,х',у',2'). 


Let ШЕТ) ШЕ ехр | wt 


ве г 
ERN oe |. 2R' 2R  2R' 
A c oc) EXD M MS 


= 
|| 


іше the energy is transmitted in finite amounts of 
EN it during those finite lengths of time the exponential 
Ш 8 по vary too much, then during a transmission the expon- 
Ва part can be considered constant although varying between 
 П501551015. So it can be taken out of the integral. The 
КЕБШЕ will be 


= _ 2Ry os . 2R' а ORE 
ж DEC - c) g*(t т rca 


EN is constant during a transmission period, it can be 
moved out of the summation and 


2R _ 2R! 


Y = [Sgg*dt] E exp -ju 1 - “С 


EN obvious that the first factor is the range ambiguity 
function, and the second, the azimuth ambiguity function. 


So, the azimuth ambiguity function is 


2R Ар 
Vaz Бо ехр -) Do ore 5 C ) 
Figure 28 indicates a real situation. From that Zısure 


ШЕИ be seen that: 
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воно му С 





_ т. x 
R e ех > Ко Е ON 
О 
EN (x Te KR + x 
О 
ШІ х << Ro k - х') << Ry: which is basically the same 


Meeeraction given to 0 in Fig. 25. 
Meche radar transmits with a PRR of fp = 1/T, then x = 


аге the positions where transmission occurs. 


which infers Vaz ӘЗ” о ES : 3) 


|| 


5 ехр - "m Tes EE S) 
О 


С N/2 | 
О Y exp -j 4т(х'/АВ ) пуТ 
c Ro -N/2 O 


exp * j 


DW Sb. the synthetic aperture length is 
SN E, у = aircraft speed. 


olx.) sin (NAAA VIT ZA RO 


ПИШЕШ infers Иду = €XP J Кос sin 4тхтут КЕ 


1122 
where exp j шо Ох”) 
Ко С 





Soma presentacion Те 
Sun (N+1)ua/sin a is the amplitude factor of interest. The 


Eun resolution is given by the 3db points separation. 
N 





158 E^ iene S db point is at 
ТХ 
= 1.4 
ARG 
1.4 A Ro 
= ! = 
and Jv 2X T 


84 





Er 


Since L 


Sl 
ө, 


NM о 


D D 
AZ Ti 2 


which is a result in agreement with that already derived. 

ншиоспегаєсіпо the synthetic aperture radar, in order to 
reduce the side lobes of VE it is sometimes useful to weight 
Bu Ereturned signals before combining them. In the case of 
Mies tocused antenna, phase compensation must also be used. 


Е № returns are processed in order to get 
L = NvT 
ШИЕ focused processing will be of the form 


LS [exp 3 Фах 


ШЕТПЕ unfocused processing of the form 
LS М 
n n 
where m the weighting function and 5, are ihe discrete 
mens. This is the basic processing technique. Various 


EENUSds, from optical, electronic, as well as З СПО ее 


used to generate the synthetic aperture radar. 
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V. RECENT DEVELOPMENTS IN TWO MAJOR AREAS 


ЕЕ DIGITAL MTI 

pire structure of an MTI processor clearly points into the 
Ex etron of a digital implementation. The increase in digital 
ЕЕ rates associated with the cost decay of present digital 
Exuponents led to the implementation of reliable, flexible, 
MnC ost digital MTI processors. 

ШЕ 29 15 а typical block diagram of one of the types of 
Коло tal implementation. I and Q channels are used in 
E cro not to lose 5 db on the average due to blind phases. 

ШОШО тис тиге is parallel to the analog MTI processor. The 

use of A/D's and D/A's, storage devices instead of delay lines, 
Bun the differences with respect to the analog processor. 
The capacity of the memory depends basically on the number of 
meee Cells and on the use of multiple cancellers. 

Mie quantization process introduces a new dimension to the 
ENS problem: that is, the quantization noise. Quantization 
norse is present along all the dynamic range of the processor, 
Ne errors at the extremes are of particular importance 
(ГС MTI processor. At the lower level, clutter cannot be 
ЕТЕУ са below one bit of quantization; at the upper level, 
abrupt clipping distorts the signal creating additional noise. 
This extra noise affects the improvement factor (I) of the 
065501. If it is assumed that signal and quantization 
noise are uncorrelated, then for a single canceller, the 
improvement factor is [14], 
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BITS RANGE GATE 


DIGITAL 
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where Em 


ЕС эргышт) AS P 
EN (о Рі | - ро (Т)р "о 3513 
DA correlation coefficient of quantization error for 


Б 
clutter alone 


> correlatıon coefficient of quantization error when 


both clutter and signal are present 


o * = noise power 


Е 
Е = signal voltage 

E. Бери power clutter 

ENS. correlation coefficient of clutter signals. 


Since ог nr = о and ре (T) ~ 0 


Ш ~ 1 


] 
1- P(T) * og7Pic 


ет тет from the equivalent analog factor only by the 


eee 
2/p. 
term o £7 Pic: 


ВЮ, or a double canceller, the improvement factor I, 


15 14 | 
E JL 
и = 4 1 2 
EE ЧОО) * « p(2T) + o 7Pic 


Meee triple and higher order cancellers, the expressions 
More complicated but the derivation is the same for all 
ENS Also іп в Peony dittLerence ronk he analog 
ОООО ШС ло js the factor А. Assuming a Uniform Error 


weri bDution 
2 En? 


O n- 
EE - 1/2)* 


Бавптаттоп voltage of the digital register 


Humber Of bits used in the quantization. 


n 
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Р с/о. is called the improvement factor limitation. 
une 50 relates the improvement factor limitation with the 
pu number of bits for different values of ek 

апу applications the fact that the MTI processor has 
blind speeds at multiples of the PRF is completely undesirable. 
iijemmethod usually used to attenuate this effect makes use of 
zc cered PRF. The scan-to-scan staggered PRF is less efficient 
but requires less complex hardware than the pulse-to-pulse 
stagger, but for some applications the first one is inadequate. 
lx taspspered PRF technique is based on the use of different 
Spacing between transmitted pulses which, at the receiver, are 
mu ly delayed (de-staggered), such that the spacings at 
HIEUDput of the delay line cancellers are all equal. The 
Me@essary use of delays makes staggered PRF perfectly matched 
Шиши Са! techniques. Since the concept is completely 
ПИ the recent trend has been focused in the development 
(шиита галоп techniques in order to improve the performance 
MN processor. The first difficulty with any optimization 
ЕЕ ЕП 15, Of course, the definition of optimum. Some try 
Пи mize the improvement factor, others to optimize the 
ЕЕЕ ТО clutter gain (SCG) within a finite number of fre- 
quency slots, or even the optimization of some defined indicator. 

Am important factor in the optimization process is the 
=ewetetical model that is chosen for the clutter power dis- 
EUH. It has been verified that a Gaussian distribution 
centered at zero doppler frequency is adequate in most cases 


Комет density function of clutter. Various standard 
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NM uscrons have been calculated for different types of 
Elutter [2]. 

Figure 31 is a block diagram of a staggered PRF processor. 
mmempulse Staggered sequence is such that the ith palse is 
memayed (1-1)T + AT, seconds respect to the first pulse. 
EN processor, after de-staggering, the pulses are weighted 
(М. "5 Ро и. Ву observation, the impulse response of 


Be processor is 
h(t) = Wy DENIS АТ.) + W, О = (Ту АТ,) ] 
eo. OE UN 6(t - [(N-1)T + 47,1) 


Euch implies a transfer function 


N 


H(w) = > M у т) 


Ша filter power response 


G(w) = H(w) H*(w) , = eonplex, conjugate. 


meen Glu) and the power spectral density of the clutter, 
E rent optimization criteria can now be devised. One factor 
is common to all techniques. The A.'S, Wi's SiGe eG кет 
variance, assuming a Gaussian distribution, are the parameters 
that will optimize the processor, given an optimization pro- 
все. In the analysis, it is always assumed a step scanned 
Eu c untenna in order that for each position, N pulses are 
ШОШ штеп. Figure 31 1s a processor for a given range cell. 
Ц го аге К large cells, the general block diagram is 
Mieaemted on Fig. 32, where W's and the summing network 1s 


presented for the К; зс. 
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iemidetietiitetulter is a high-pass filter with rejection 
Шо лето frequency to tlle maximum frequency of the clutter. 
Mepectionemay also occur after the maximum desired doppler 
Frequency but it is not absolutely necessary. 

Шеп а Detter understanding of staggered PRF some general 
ниееро5 should be explored. If a constant PRF is used, not 
only the dc component of the doppler spectrum but the whole 
EDU cr spectrum is translated to multiples of the PRF. It 
Exon [l5] that with staggered PRF the clutter spectrum 
Ems translated to multiples of the average PRF, but the 
spectrum of the signal is dispersed into N (number of A.'S) 
EN cte frequency lines. So, if a target produces a doppler 
mea uency that is a multiple of the basic PRF, only one part 
ШІ (Пе whole signal is distorted by the clutter spectrum 
Dhere the remaining M-1 parts of the signal may possess enough 
meto be detected. An MTI filter which notches in a region 
Bud dc as well as in multiples of PRF can eliminate the 
ИЙЕ СТ, but if staggered PRF is used only one part in N of 
ИО ОТОП 5 1б65$С, thus increasing the probability of detec- 
EO targets with potential blind velocities. 

№0 16| an optimization process, based on a performance 
ЕР defined as 

B 


а ағ 


51 [s2(£)]" 


EE veloped. B4 and В, are the lower and upper bounds of the 
LN uWrv region, S?^(f) is the signal to clutter gain and M is 


Emeter that reflects the emphasis the criteria puts on 


э? 





spectral amplitude variations. As M increases, smaller 
ше ир гЕПасс have more influence in the final result 
Bi спе integral. 

Oia derstanding of the significance of P, it should 
be noted that good clutter attenuation as well as a high im- 
provement factor correspond to a good performance index. 
ЖР ЕПС the optimization procedure described in [16] and using 
practical results (figures 32, 55, 34) for different 
values of N and (Ба - By) were obtained. It can be seen that 
increases, the minimum value for the signal-to-clutter 
gain in the velocity region increases as well as the improve- 
pent factor I or acp Dat with a fixed N and 
Пеј ослту region increased on the upper bound, the minimum 
memes Of SCR(f) in the velocity region decreases (Fig. 33, 34). 

ша sala earlier, the clutter variance ø is also an 
eme tant parameter in the optimization process. Fig. 35 
дише о With the improvement factor for the case of N = 5 
and В, = 18 Ву. О ОЕ ИВЕ е interpulse periods 
Rt A; were optimized for the value of o = 0.02 Ву and the 
пе слеп 5 М. were Optimized in parallel with the variation 
Ma the horizontal axis. On the dotted line, both coef- 
ficients and interpulses were optimized for o = 0.02 Ву, апа 
NE r1atrons of I with co were recorded. This situation 
occurs when the real clutter parameters are not those used 


Mato design of the filter. It is apparent that there are 


not too many differences between the two cases. 
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SCR gain as a function of Doppler frequency; 
N =4,B,,=108,,0=0.028, . Optimum intervals: 1.100, 


1.094, 1.000. Optimum coefficients: 1.000, —3.003, 
3.155, —1.152. 


Fig. 32 
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SCR gain as a function of Doppler frequency; 
N 5, 8,, = 108,. 2 = 0.028, .„Оритит имегма!5: 1.111, 
1.000, 1.091, 1.058. Optimum coefficients: 1,000, 
—4.258, 6.242, —4.031, 1.048. 


Fig. 33 
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SCR gain as a function of Doppler frequency; 
N=5,B,= 408,,0 =0.028, . Optimum intervals: 1.007; 
1.021, 1.000, 1.299. Оритит coefficients: 1.000, 
—3.784, 5.432, —3.273, 0.626. 
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Improvement factor as a function of clutter spectral width 
oi — 5, B; = 108, . Coefficients for dotted curve (optimum for g = 
0.028,): 1.000, 4.258, 6.242, -4.031, 1.048. Intervals for both 
solid end dotted curves (optimum for y = 0.028, ): 1.111, 1.000, 
moo}. 1.058 (С 1. Рао. 9). 
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Another criterion to optimize the staggered PRF processor 
ШОС oped in [17]. Instead of maximizing the average sig- 
e pain Over the entire velocity region, the 
ШООК Осиет gain is maximized in small Af intervals 
ев the velocity region. The number of intervals is, of 
не гри тесту related to the complexity of the processor. 
шп the optimization algorithm developed in [19], a final 
ЖИБЕК diapram for the processor is reached (Fig. 36, 37). In 
ео спе input signal after being split in I and Q channels 
Ше мелеһісі with the optimum set of weights A and B for each 
Пише. The results are then combined in order to get an out- 
MO positive doppler as well as an output of negative doppler. 
EE represents the implementation of the filter that gener- 
EN us sets ot A and B weights for each of the M frequency 
sub-intervals in a channel. The Т.'5 ape the interpulse delays 
апа the 

"uu |... UN c 1 

і 0<j<M-1 

ehe coefficients determined by the optimization process. 
should be noted that in this optimization process the 
EM rpulse periods are not optimized: only the weights are 
me@emnized. Fig. 38 relates the signal-to-clutter gain aver- 
aged over the entire velocity region with the signal-to-gain 
maximized in each of eight frequency sub-intervals. Both 
N = 5 and N = 10 are presented. A substantial improvement 
due to optimization is evident. Another optimization tech- 


ме 15 developed in [18]. 
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Mia cd processors are called Constrained Improvement 
INP roces sors (СТР). Given a specified improvement factor, 
the number of pulses to be processed and the PRF stagger 
Sequence, the mean square deviation from a constant response 
MNG velocity region iS minimized. Two main approaches can 
ШО ento the problem. In the first, the PRF stagger 
Sequence js fixed and the weighting coefficients are chosen 
EN xooprtimrze the processor; in the second, the weight- 
Mmg coefficients are fixed while the stagger sequence is chosen 


order to optimize the processor. Defining 


maximum PRF 


Hh 

| 
"Tj cer 
| И 


unstaggered PRF 


Цео 59, 40 and 41 represent the optimum filter response 
[maimed опе ри] ѕе return, £' = 8, go = 0.01 and [I = 30 db, 
Meme respectively a linear PRF stagger with +20% interpulse 
Variation, and a sinosoidal PRF stagger with +10% (Fig. 40) 
me 90s (Fig. 41) interpulse variation. As it can be seen, 
EN ne variation increases the response becomes more uniform. 
Eupsrunpg figures 42, 39 and 43, in which the number of pulses 
are respectively 3, 4 and 6 for a common interpulse variation 
Of 20$, it can also be seen that the response improves with 
спе number of pulses processed. 

A comparison between optimization procedures, different 
from a direct comparison between responses, can be devised 
Meme the fraction of frequencies for which a response is 
less than some specified value as a comparison parameter. 


meee 0 db scale corresponds to the mean value of the 
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Frequency response of four-pulse CIP using sinusoidal PRF 
stagger with £ 10 percent interpulse period variation, 
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frequencies of interest, that means that 50$ of the frequen- 
cies are above as well as below the 0 db level. This leads 
to a cumulative distribution function in which instead of a 
Нол о: trequencies, the term probability is used. Fig. 
Esca cumulative distribution of the responses of a four 
ШШ ОСОО СТР апа a four pulse processor using the Ref. [16] 
пита блоп стабетја. This comparison criteria gives a 
member performance for the CIP processor since its curve is 
Moser to being ihe ideal step response at 0 db than the other 
one, but on the average they are similar. 

| Another песни шоа MITT processor is called the matrix 
fee in this case there is no parallel between the analog and 
Па тах MTI processor. This processor (Fig. 45), in addi- 
EN ue tremendous flexibility in the modification of the 
Shape of the response curve, has associated with it a very 
Ке Clutter locking mechanism that shifts the response 
En cwhenever there is an average velocity associated with 
NESCIT. Іс simplify the explanation of the processor in Fig. 
45, опју two levels of quantization are used. Basically the 
 ЕЕС550Т7 15 a coherent System that compares the phase differ- 
ence between two successive returns and weights that difference 
wae according to the shape it is intended to give to the 
response curve. Since dy = 044%, by weighting dd in effect 
 борріст 15 being weighted. 


ЕЕЕ ШӘП Ос а Signal return is defined by 


о = the > 


101 





CIP 


e wv PRINSER 


38700 -20.00 -10.00 0.00 10.00 20.00 
КЕ5РОМ5Е(08) 


102 





ino 
1и 





SAYOMLIN 
9NILHOI3A 


U3LNNOD 


Ep) 
< 
oD 
r1 
E, 








83110112 v 3o 1l 
01 
пес 
E Y 
43191935 © 
XIYULVW u 
NOISSIUINOI i 
^ 1и : 1 
En 1182 | | 
~ малбозн | HS 13 
07 
eee eee! ee Ce TO EM. 
| 
ual 
| Rn v/a : 434743 b] 024 | 
| | 
| | 
| 1105312 Я201 8311012 | 


ТҮ М215 
азил 


103% 





and the amplitude by (I^ + 0?) | О сите саге of 

INIM по се only once quantization level for I and Q, 

that means that the only possible phases for I are 0° and 
Мапа for Q, 90° and 270°, assuming a signal to noise ratio 
much greater than 1. When the average phase difference is 0°, 
Што 1е5 а <та опагу target, but if doy = 180° that indicates 
mun velocity target. Since each vector return can be 
Теа in all of the four positions ЧІ = jQ, there are 4x 4 
поје Combinations for two consecutive returns. So, what 
ши оотруе55топ matrix block in Fig. 45 does is to generate 
ита 15 proportional to the phase difference output dor 

E cen two consecutive vector returns, which in this partic- 
Menase Can take the values 0°, 90°, 180° and 270°. Each 
ene four outputs мәсе а Gomme sponding weighting factor W 
ЖИЕ ПЕ result is summed to generate the MTI output. Fig. 46 
E uu ces the four types with the vector combinations that 

ШО ласа them, as well as the weighting factors used in order 
И оотохтасе an MII response curve. The final response curve 
Semresponds to a statistical average of d$, when uniformly 
EXON phases are introduced at the input. In order to have 
Pees! proportional to the clutter velocity, it is only 
mecessary to subtract the outputs of the region (d) (Fig. 46) 
ШЕШ (ПоО5вес oft region (c) and divide by the number of range 
EN that is, for two quantization levels and a uniformly 


ШЕ i buted phase input 
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So, шп а large clutter environment if N is large, do, can 
smoothly adjust the local oscillator to place the MTI null 
ЯС (пе mean value of the clutter doppler. 

Поле се mew ann Quantization bit circuit. For each 
channel there are 2” poe De paguant ization levels. If only 
Smemouddrant Of possible phase differences is considered, it 
can be seen that there are (Fig. 47) 

ERU T a, 


2 2 


к Шрагсеі results from the fact 


possible phases. The -2° 
that of all the diagonal combinations (45° zone) only one 
combination can be counted. As an example, let's use two 


 шГгігатсіоп bits. That means that there are 


uem n-1 


4[2 12 


~ 
+ 

ER 

md 
i 


п= 2 


possible phases in all four quadrants, which also implies 12 
possible phase differences. 

ШІРПЕЕ 46 (а) represents the response and figure 48 (b) 
Excunatrix output (phase) for each 12 x 12 possible com- 
Eu con Of two consecutive vector returns. Since there is 
EM ubility in modifying the weighting coefficients, it is 
possible to adapt the filter as a function of the type of 
Expected clutter spectrum. Again, the mean clutter velocity 
can be determined by subtracting the outputs of the (b) region 
(Fig. 48) from those of the (m) regions and applying the 


meee tO the VCO associated with the local oscillator. 
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Total number of combinations per quadrant 
Salta tle diagonal combinations only as 
Se 2-1. 21-1. 21-1 К 

Total number of combinations for the four 
quadrants - Eo ml 4n-l , B 
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P DIGITAL TSAR 

The synthetic aperture radar (SAR) is a typical example 
расла processing where digital techniques are being sub- 
De edi or the former optical processors. With the present 
й знпотосу mini-computers can compensate for the aircraft 
movement, and, as with a real antenna, the synthetic antenna 
ШЕП БС steered and even scanned. The basic problem with SAR 
НС сссок5 15 the need for high storage capabilities and very 
fast data rates. The film as an optical storage device became 
most impractical since a real time display is impossible. 
Meuse Of storage tubes is also inadequate due to low effi- 
Еа псу and poor dynamic range and stability. With the present 
high speed A/D converters (100 MHZ), low cost, and high speed 
Exact Gigital storage devices, it becomes feasible for real 
meme SAR processors. 

Ке У shows a time Overview of the late achievements 
SAR techniques. 

Figure 49 is a block diagram of a digital synthetic 
E rure radar. Theoretically, the processor must, for each 
те gate, perform the integral 


Ос (Ст аът 
АТ 


Were S(t) is the signal return and r(t) the correlator 
erence function. 
r(t) = Ap(t) exp -j óp(t) 
Ap (8) IS ЕИО LON to control the side lobes 
othe synthetic antenna pattern. 
dp (t) = phase reference that tracks the phase of s(t). 
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1957 


Mid 1960's 


Late 1900" 5 


волани 197075 


Е Development 


Carl-Wiley postulates doppler beam 
ӘНШІ со псе рт. 


Dyer iy Ot eines demonstrates 
ШЕШ ет репе concept. 


Ici Кн сокири ло optical 
cone ao ale proa 


ОСИКОЮ К Се sonk correlation 
demonstrated in non real time. 


ПП КЕЛЕСОО ТОП ТС SAR correlation 
demonstrated in non real time. 


ПОСНА о на AR demonstrated 
CASO CON pens at 101. 
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Зіпсе фу О) follows the phase of s(t) but with an unknown 
@umreerence, the Signals are processed in an I {in phase) and 
EN Out Of phase) Channel in order to prevent signal losses. 
NcEdyprtal processor operates in the following way: The radar 
returns pass an A/D converter that samples the signal at a 
moe, at least higher than the Nyquist rate, and separates 
ШШЕ гейптүпс in range bins which are approximately equal to 
Пе апре resolution. The number of bits used per range bin 
Кеа function of the desired dynamic range of the processor. 
The digital data passes through a buffer and prefilter which 
ПИ аге the A/D) rate to the rate of the correlator. The 
mee memory stores the data corresponding to an integration 
time, AT. The data is then correlated with the reference 
peveratOr output in order to form an azimuth line with a 
length corresponding to the number of range gates. After 
Mo lation, new data enters the memory and the process re- 
uus itself to produce a new azimuth line. 

ШЕГІ 3) 1! сап be Seen that the bulk memory in bits 
ШЕГЕ с сапа! to the number of cells times the average number 


Seats per cell, so 


BM = bulk memory = ¿Ko М Nn 
E : _ AR 
Np = number of range gates = = 
Y 
N = number of azimuth data lines = AT = ШУТ 


Т r 
uu uctor of 2 is because of the use oi an 1 and Q channel. 
Et. for example, the desired range resolution is 50 m, 


Bat corresponds to a bandwidth 
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BW - >= ВТО м сес _ 
а 2 о и 


O - 

which implies an A/D sampling rate per channel of 5.0 MHZ. 

On the other hand, for a 1 KHZ PRE the unambiguous range is 

80 mi, which for a typical AR = 10 mi of mapping, means that 
ПИ легла | ап equivalent 70 mi excess time to process that 
Шта. This results that in fact there is no need for the 
@erreiator to work at such high data rates. Using a buffer 
ШЕН ассертв the data at a very high data rate during a short 
mod of time but delivers it at a lower rate to the corre- 
ШОО ОПГ 25 possible to reduce the correlator rates by many 
Orders of magnitude. In fact, analyzing the function of the 
mememetator, it can be seen that it must in AT seconds produce 
Amber Of Outputs equal to the total number of cells, that 
NS. (L/R,)N,, where КЕ = СА 15 the azimuth compression rate. 
Since the real antenna illuminates (L/r,)N, cells and is Teady 


to receive new data after T seconds, the correlator rate (CR) 


is 
св = (/Та)№ = САМ, # 
і і _— Rm 
T 

defining 

Мы - Ка Mabe tor: doppler filters 

De. Не ОВ sample factor = 1 
Then 


СК = Ny Np f. 


Yt Np = 500 and t. = 1 KHZ for one doppler channel CR = 


Benz < A/D sampling rate. 
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шев вв Қапа, 11 the Spacing corresponding to V/T 
iS much finer than the range resolution Гү, that means excessive 
data is being used. If the PRF cannot be lowered due to power 
BENEEdODDpyer considerations, then a prefilter can be used to 
effectively reduce the PRÍ = Eno by a factor f/f, where EN 
SN ne prefrlter sampling rate. The prefilter will reduce 
NU Scorrelator rate by a factor БН нет те итге т е ше 


Mon 15 in order to have V/T of the order of Po then 


N 
I ~ EL. EE ке > ис а nEw Correlator rate. 
5 m ra Г Т 
= E 3 2 
ER Sa Ko МЕ No f. 
ie ou оле арс Factor 
ko = synthetic array weighting 
MU Sesto factors, the bulk memory and the correlator rate, 
determine the type of design for a SAR processor. The main 
ШИССІІУСс 15 to reduce both CR and BM. Today's technology is 


ЕТЕ ос this problem in two ways: devising correlator algor- 
вепт5 which will reduce the bulk memory and arithmetic, at low 
E Сп fewer power consuming elements. The first approach 
meebeing made by parallel and series combinations of correlator 
@mammels and prefilters, as well as with the use of FFT algor- 
meee Another method is the use of pulse compression tech- 
niques but processed only after the SAR processor. This will 
T educe Ка which will reduce BM. The hardware improvements 
mun CCD memories and LSI at much lower costs. 

Pr ла said before, with today's digital techniques 


various mapping modes are possible as well as motion compensation. 
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Pigure 51 shows the four types of mapping used today. The 
МИНИ осе (Ра. о: Б) 15 опју а variation in angle of the 
side-looking SAR. The doppler beam sharpened mapping is a 
Pelerepresentation with an equivalent antenna of very high 
azimuth resolution. ЩЕ О КЕБ а mapping where a high 
resolution snapshot map is generated. Motion Compensation 
NUM dos various functions: Clutter tracking, focusing accel- 
Beton compensation and real antenna stabilization. Fig. 52 
EN NSblock diagram of a general SAR processor. The blocks 
ме разаса ју the same as those of Fig. 49 except for the 
motion compensation blocks which are now included. 

To achieve lower data rates, various types of algorithms 
are presently used. Lower rates will be achieved at the 
pense Of more complex hardware. The general concept that 
ИЕР ТО dali of them is the reduction of rates by the in- 
Gmease in the number of doppler filters. 

(то о ista block diagram of the multi-channel prefilter 
Processor approach [19] where m is the number of channels. 
Ри отет] бег rate becomes 


Вр = m Np i 


and since the number of filters Ne is reduced byl/m, the total 


arithmetic rate (TAR) is 





LM DR. вт. 2 
a, = No f. 
а,= Ко. К, Мр Np/AT 

d TAR Фо _ 
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M MESEchgt o the vatue of m that minimizes TAR is 


1 
п = [%2} = [К 


р 2 1/ 
К. Np/£, AT] 2 


OS 


ad the minimum TAR is 


TAR 


min 


С С, 


у 
= 2 T 
2(K oe K. f./AT) М.М, 


OS 


Mies two-stage correlator (Fig. 54) is another algorithm [19]. 
The idea is again frequency domain division but only with one 
И | со lie eguivalent correlator rate will be the sum of 


meme two correlator rates. So 


ТАВ = (Ко. K. No N¿/0T) M 


i ћ 
+ Kos as ES ШІ 


where M is the number of first-stage correlation channels. 
As in the previous case, there is a value of M = (N,) ^ that 
ШЕШЕНІ 765 ТАП. With the FFT algorithm, it is possible to 
luce even more the total rates, Table V is a summary of 
Emezealculated values [I9] of bulk memory as well as total 
EE metic rates for the SAR signal processing algorithms. 
mother approach, completely different from those already 
ши оепоед to solve the problem of high data rates on the 
@ammeclator (Fig. 49) of the SAR processor, is through the use 
GE parallel processing using associative memory [20]. Fig. 55 
represents a block diagram of the associative memory and the 


Meme ted Input/Output registers. Іп the associative memory, 


T19 





ЗЕТ 


2801035 10189102914 о) әш? ола 156 ay) PRU 29301015 101018910! 20] 97e surpirrogpe doy) ISM ә] "10121309401 әу ЧІ 10 1012191 
-102914 51 280105 оці гоціоцм цодп 5итриәдәр ‘ospe "рџе зшу)нозје ou Suoure Apsi[s A174 [rt Py s31Q OPn3riduse Jo soquinu oy :2:0 


> 


DEA IN А Y зе 144 518 Јәиуоі 


Uno Рн /| (үйме кд АЛЗА УЗО уус | оу одејо-ома ща ән 











МА Nz UV AAA nz 103c91100 snid s123[21d op£tmy 
AS з) = vn о APT AP xt 1012:31103 51,4 1911Ј234 
UT V WA xz 107243110) 

ILY опошцију [eor Alowaly NIN шутшощу 





ЕЕС m т. 


120 








eS зі 
МАСС БЕСІСТЕБ | 


noe VE 


ADDRESS > 
| MEMORY E > 
REGISTER 
REGISTER 





Bee 1 


OUTPUT 


о 55 





Na dres sing, multi read/write, multi logical and arith- 
metic operations are performed. The search/write and the 
DExcesysterscare used to generate proper codes for multi 

Ep ssane. The operations are done in a bit serial basis but 
on all words in parallel at the same time. Ellis [20] showed 
mut for an X band radar with range and azimuth resolutions 
of 10 meters, 1000 range cells, maximum range of 150 km and 
Maximum aircraft velocity of 200 m/sec, a 50 msec correlation 
Ens necessary: this implies a 40 nsec multiplication time. 
Mees associative memory, it is possible [20] to perform the 
шта јоп5 OL a correlation period in 31.48 msec. The advan- 
mus of hardware parallel operations are thus obvious over 


the conventional process when speed is an important factor. 
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КИНО ОДА GT ONS. TRENDS 


It was shown that due to improvements in integrated 
ШОО ту better macroprocessors, faster DFT/FFT algorithms 
and lower cost of digital logic, most analog processors are 
Deine replaced by digital processors. Very high data rates 
пе to be the major problem of digital processors, but 
ЦЕН (ойау |5 techniques of parallel processing, general pur- 
pose digital radar signal processors are already in use. 
Mespite the quantization noise inherent to digital processors, 
ШОО reliability, extreme flexibility and low cost of digital 
MmecessOrs Still give them a tremendous advantage over the 
mermer analog processors. 

EE cture demand for more reliable and sophisticated 
medar systems will be a function of cost and military neces- 
ЕТТІ the cost of digital logic trends lower as is its 
село trend, and military demands continues high, digital 
meessors will play an increasingly important role in modern 


Ш аа system's implementation. 
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